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Abstract. We describe several results concerning global quantum quenches from 
states with short-range correlations to quantum critical points whose low-energy 
properties are described by a 1+1-dimensional conformal field theory (CFT), extending 
the work of Calabrese and Cardy (2006): (a) for the special class of initial states 
discussed in that paper we show that, once a finite region falls inside the horizon, its 
reduced density matrix is exponentially close in L 2 norm to that of a thermal Gibbs 
state; (b) small deformations of this initial state in general lead to a (non-Abelian) 
generalized Gibbs distribution (GGE) with, however, the possibility of parafermionic 
conserved charges; (c) small deformations of the CFT, corresponding to curvature of 
the dispersion relation and (non-integrable) left-right scattering, lead to a dependence 
of the speed of propagation on the initial state, as well as diffusive broadening of the 
horizon. 


Quantum Quenches to a Critical Point in One Dimension 2 

1. Introduction 

In a global quantum quench, a quantum system extended in space is initially prepared 
in a pure state | "0o) ? which is usually taken to be the ground state of hamiltonian H 0 
with short-range interactions. It is then evolved unitarily with a hamiltonian H , which 
usually differs from H 0 by the quench, or sudden change, of some parameter. Such a 
protocol is an idealization of real experiments which may be carried out, for example, 
on ultra-cold atoms. However, for the purposes of this paper, the important features of 
|-0o) are that it is translationally invariant in the thermodynamic limit and that it has 
short-ranged correlations and area-law entanglement entropy. Furthermore we consider 
only the case when H is gapless, corresponding to a quantum critical point, and with a 
linear low-energy dispersion relation oj ~ v\k\. In that case the low-energy spectrum is 
described by a conformal field theory (CFT) with hamiltonian Hcft- 

A great deal is known about the properties of CFTs, especially in 1+1 dimensions. 
In 2006 Calabrese and this author [Ij (CC) used this knowledge to derive properties of 
a quench to a quantum critical point whose entire spectrum is described exactly by a 
CFT, that is assuming that H = Hcft- They also made a particular ansatz for the 
initial state I'i/’o)- Their main results were: 

• 1-point functions of non-conserved local observables in general decay exponentially 
to their values in the ground state of H ; 

• the ‘horizon effect’: two- and higher point functions of local observables whose 
arguments lie in a finite interval A of length i become stationary after a time 
t ~ £/2v when the left- and rightmost points first fall into a forward light cone 
whose apex lies on t — 0; 

• thermalization: after that time they are equal, up to exponentially small 
corrections, to the values they would have at the finite temperature /3 _1 in the 
CFT corresponding to the mean energy E 0 = { , iPo\Hcft\'^o)■ This is equivalent to 
the statement that reduced density matrix of the interval pt oc Tr^ eT ^ HcFT ; 

• although the whole system remains in a pure state, the entanglement entropy of the 
interval with the rest of the system increases linearly with t up to the equilibration 
time, after which it becomes stationary and oc t\ the extensive part is equal to the 
thermal Gibbs entropy at temperature /3 _1 [12]. 

• the horizon itself has width 0((3), and it is only inside this relatively narrow region 
that the details of the particular CFT become relevant. 

Although these results depend on the specific assumptions discussed above, as 
well as detailed analyticity properties of the CFT in complex space-time, they may be 
physically interpreted rather simply in terms of the emission of left- and right-moving 
pairs of quasiparticles from the initial state, which are correlated and entangled over 
only lengths 0(/3), and which otherwise move semi-classically. It turns out that, suitably 
adapted to account for the dispersion relation, this quasiparticle picture extends to other 
hamiltonians, including lattice models and to non-critical quenches, at least as long as 
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the hamiltonian H is integrable so that even the high-energy states have a quasiparticle 
description. The horizon effect is supposed to hold more generally in most systems as 
a consequence of the Lieb-Robinson bound [3], and has been observed qualitatively in 
experiments [Tj. The general picture of thermalization has also been discussed within 
the AdS/CFT correspondence, when it is supposed to be a consequence of the formation 
of a black hole [5], 

However this leaves several unanswered questions, some of which we address in this 
paper: 


(i) can we quantify how close pi is to a thermal density matrix? 

(ii) in an integrable model, there are an infinite number of conserved quantities which 
commute with the hamiltonian H. In that case, it has been suggested that 
stationary state should be described by a generalized Gibbs ensemble (GGE) 
rather than a simple Gibbs distribution [6], This has been the subject of intense 
investigation in several solvable models [7]. How does this appear within a CFT, 
which (at least in the rational case) is the most integrable model of all? 

(iii) what happens if we relax the particular ansatz that CC made for the initial state? 

(iv) what happens if we add irrelevant interactions to Hcft which give only non-leading 
corrections to the low-energy behaviour? How do they affect the quench dynamics? 

(v) what happens in a finite system? Is it possible for the initial state to fully or 
partially revive? 

In what follows we shall show that, under the assumptions of CC, the answer to (i) 
is that the overlap between pe and the reduced density matrix pp in a thermal ensemble 
is exponentially close to unity, that is 

i- K const e -2»A((-(/2„)/|) (1) 

((Trpf)(Trpg)) / 

where A is a universal exponent depending on the initial state. (The left hand side is 
of course always > 0 by the Schwarz inequality.) 

The resolution of (ii) and (iii) is as follows: the ansatz of CC amounts to the 
assumption 

|Vh) tx e~^ /4)HcFT \B ), 

where | B) is a conformal boundary state. (The factor /3/4 is chosen so that the 
expectation of Hcft hi this state is the same as that at finite temperature /3 .) We 

remind the reader that each bulk CFT is supposed to have a particular allowed set 
of such boundary states, each of which corresponds to a fixed point of the boundary 
renormalization group (RG), each with their own basins of attraction. For the minimal 
models, there is a finite number of these states, and their basins of attraction are 
supposed to contain the ground states of each possible non-critical hamiltonian Hq. Thus 
this ground state should be representable in terms of all possible irrelevant operators 
acting on | B): 

IV’o) oc J^e-ft/^15), 

k 


( 2 ) 
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suitably regularized. One of these irrelevant operators is always the component T tt of 
the stress tensor, whose space integral is the hamiltonian Hcft■ Thus the CC ansatz is 
equivalent to assuming that this is the only term in the above sum. T u may be written 
as a sum of holomorphic and antiholomorphic operators T and T, each of which is 
conserved and whose space integral is a conserved charge. This is equally true of all the 
other irrelevant operators which are its descendants, which can be written as powers of 
T and derivatives thereof, plus their antiholomorphic partners. We shall argue that in 
these cases the quench dynamics from such a state with a hamiltonian H = Hcft leads 
to a stationary reduced density matrix for a finite interval which has the GGE form 

p e oc Tr^ JJ e~ h cos(?rAfc/2) f(*k+*k)dx , (3) 

k 

where now and 4+ are holomorphic and antiholomorphic bulk operators, so that their 
space integrals are conserved charges. Note, however, that the charges f (4+ + §k)dx 
above do not necessarily commute among themselves, even though they commute with 
the hamiltonian Hcft■ The conventional GGE includes only a commuting sub-algebra 
of these charges. This is motivated by the idea that these should form a complete 
set of macroscopic commuting observables which should characterize any macrostate. 
However, this does not appear to be the case for a 1+1-dimensional CFT: because of 
the exactly linear dispersion relation, there is a massive degeneracy of states, and the 
expectation values of the charges in the commuting sub-algebra (which were identified 
in [8]) are not sufficient to characterize the states. A similar phenomenon has been 
pointed out by Sotiriadis [9] in the context of a quench in a 1+1-dimensional massless 
free boson theory from a non-gaussian initial state, although in this case the commuting 
conserved charges are the mode occupation numbers which are non-local. 

However, this still does not account for all possible irrelevant boundary operators 
in (J2j). We further argue that, at least in rational CFTs, for every boundary operator 
4+ there is a pair of holomorphic and antiholomorphic operators (4+, 4+) with the same 
overall scaling dimension A*.. These also lead to conserved charges and a generalization 
of Q. However, since A^ is in general non-integer, these currents are parafermionic and 
the corresponding charges are non-local, acting, for example, to change the boundary 
conditions. Although such charges are not customarily included in the GGE, we shall 
argue that they make perfect sense, at least perturbatively in the /3k, and that they lead 
to physically sensible results. 

Question (iv) seems to be difficult to address in general. A perturbative approach 
would require knowing the correlators of all possible irrelevant operators after the 
quench. However, in the case of operators (T 2 , T , TT ) involving the square of the stress 
tensor, which are present in all critical theories with a UV regulator, we can express the 
effect in terms of a coupling to a random metric with short-range correlations. These 
terms may be thought of as introducing soft scattering of the original CFT quasiparticlcs. 
The first two terms break the low-energy Lorentz invariance and correspond to curvature 
of the dispersion relation. The TT term corresponds to left-right scattering and is the 
simplest non-integrable perturbation of the CFT. The geometrical interpretation of these 
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terms leads to two main effects: a modification of the propagation velocity v dependent 
on the effective temperature /3 -1 of the initial state, and a spreading of the horizon 
width as (t/3) 1//2 . Although strictly valid only at weak coupling, it may be argued that 
this general picture should persist to higher orders and for other irrelevant operators 
which are powers and derivatives of the stress tensor. Dependence of the propagation 
speed on the initial state has been observed in several studies of lattice models [10] . 

(v) has previously been studied in Ref. HU- It can be addressed in terms of the 
return amplitude 

F(t) = \(^o\e~ lHt \^ 0 }\ . 

In a periodic system of length L (and also in an open one with conformal boundary 
conditions) in pure CFT the eigenvalues of H are of the form (27 t/L)(A + integer), and 
therefore in a rational CFT, when the scaling dimensions A take only a finite number of 
rational values N/M, there will always be complete revivals with T = 1 at times such 
that 2 vt = ML (vt for an open system). However, whether there are partial revivals 
with T = 0(1) at earlier times depends on details of the CFT. For the CC initial state it 
turns out that T(h) is given by the continuation of the partition function of the CFT on 
an annulus, or rectangle, continued to complex values of the modulus, or aspect ratio. 
For rational CFTs this leads to partial revivals whenever 2 vt/L is an integer, and in 
fact a remarkably rich structure near every rational value of vt/L , related to properties 
of the CFT characters under the modular group. If a weak perturbation such as TT is 
added to the hamiltonian the effect is to broaden the revivals with a width oc t 1 / 2 . 

The layout of this paper is as follows. In Sec. [2] we recall the arguments of CC, in 
particular rephrasing the basic assumption about the initial state made in those papers, 
then recall the main results, their physical interpretation, and the technical properties 
of CFT that allow these arguments to go through. Once this is set up the statement 
about the convergence of the reduced density matrix to the thermal distribution pp 
follows fairly straightforwardly. 

Sec. [3] addresses deviations from the CC initial state, and whether they lead to a 
GGE. We first restrict to the case of irrelevant operators related to the stress tensor 
and derivatives and powers thereof. We then study the properties of this GGE in 
CFT in general, and show, for example, that it leads to observable-dependent effective 
temperatures. Since these arguments require the use of possibly delicate analytic 
continuations within the path integral, it is reassuring that the results can also be 
obtained by a straightforward resummation of perturbation theory. This also leads to 
a simple physical interpretation of the results of this section. 

We then discuss the slightly more tricky question of more general perturbations 
of the CC state with non-integer dimensions. We argue that their effects may only 
be described within a GGE if semi-local charges which are integrals of parafermionic 
currents are included. We give an illustrative example of the transverse held Ising model 
quenched to the critical point from a short-range disordered state with a small non-zero 
longitudinal magnetization. 

Finally in Sec. [4] we consider irrelevant perturbations of the evolving hamiltonian 
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Hcfti specifically those second order in the stress tensor. We argue that these are 
equivalent to coupling to a random metric, and, to lowest order are able to deduce the 
results described earlier. We then argue that the qualitative picture they give remains 
valid at higher orders and for more general perturbations. 

2. Convergence to a thermal ensemble. 

2.1. Correlation functions. 

In this section we revisit the arguments of Calabrese and Cardy m pc] , reformulating 
them in a way that makes the assumptions more transparent, and then argue that 
the reduced density matrix pi of a finite interval, after falling within the horizon, is 
exponentially close to that of a thermal ensemble. 

We consider the time evolution of a one-dimensional quantum system from an 
initial state |^ Q ), assumed to be translationally invariant, with short-range correlations 
and entanglement. The system evolves unitarily for times t > 0 with the hamiltonian 
H , which we take to be that of a 1+1-dimensional CFT. Thus the state at time t in 
the Schrodinger picture is e^ ltH \f)o) with H = Hqft- We are primarily interested in 
correlation functions of local operators dq {x 3 ), 1 < j < n, possibly evaluated at different 
times, which will therefore be of the form 

(H J tnH ® n {xn)e- i{ - tn - tn -' )H ... Mx2)e-< t ^ H Mxi)e~ it ^o) , 

where t n > t n - 1 > ... > t\ > 0. Since CFT is most easily formulated in euclidean space, 
it is simpler to start from the imaginary time evolution 

(fj 0 \e- TfH ^n(x n )e^ Tn - Tn -^ H ... ^ 2 (x 2 )e-^- Tl)H ^ l (x 1 )e- TlH \fj 0 ). 

This can be considered as a correlation function in the Heisenberg picture 


{*&n( x ni T n) ■ ■ • $2 { x 2, T 2 )$i (xi, Ti)) , 

in a euclidean slab geometry 0 < r < r to t = T n + Tf with boundary conditions on each 
edge corresponding to the state |^> 0 ). If should then be continued to t 3 — > itj and 
Tf —> —it n . The problem with this is that the total width Ttot of the slab is then zero. 
In CC this was avoided by appealing to the theory of boundary critical phenomena: 
the actual boundary conditions at r = 0 and r tot are replaced by ‘idealized’ boundary 
conditions at r = —r 0 , r to t + To, where To is the so-called extrapolation length. In the 
theory of boundary critical behaviour, this is justified on the basis of the renormalization 
group (RG): the idealized boundary condition corresponds to a fixed point of the RG, 
and To measures the deviation of the actual state from this. For the example of a 
transverse Ising model initially in the disordered phase, this ideal state would be an 
unentangled product state. Once we take a finite To, it then makes sense to take 
Rot —> 0. We must then compute correlation functions in a slab of width 2r 0 with 
ideal boundary conditions. However this argument is based on the assumption that the 
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long-time behaviour after the quench should be insensitive to the details of the initial 
state as long as it has only short-range correlations. As we shall argue, at least for the 
case of evolution with a CFT hamiltonian, this is not in fact the case. 

Let us therefore rephrase the prescription of CC in a way that its assumptions are 
more clear and so may be generalized. It is equivalent to assuming that the initial state 
has the form 


IV’o) oc e~ ToH \B) , (4) 

where | B) is the (conformally invariant) state corresponding to the idealized boundary 
condition, corresponding to a fixed point of the boundary RG. In CFT, such fixed points 
correspond to conformally invariant boundary conditions, for which the component 
T tx oc T — T vanishes 

(T(x) -T(x))\B) = 0. 


For such boundary conditions, correlation functions in the slab are simply related 
to those in any other simply connected region obtained by a conformal mapping, in 
particular the upper half plane HI. 

Note that the expectation value of H in this state is given by the free energy per 
unit width of a euclidean strip of width 2r 0 

ncL 


{i} 0 \H\i/j Q ) = 


(5) 


24(2r 0 ) 2 ’ 

where L is the total length, c is the central charge, and the expectation value in the 
ground state has been normalized to zero. This may be compared with the mean energy 
at finite temperature (3~ - 1 pass] 


Tr He-? 11 _ ncL 
Tr e~h H 6(3 2 


Since H is conserved we see, anticipating thermalization, that To = (3/4 and we adopt 
this parametrization from now on. 

Altough from this point of view, the prescription (]4]) seems rather artificial, as we 
argued in the introduction it may be viewed as including just one, in fact often the 
most important, of the possible irrelevant operators in the sum (l2lh For the purposes of 
the remainder of this chapter we shall assume this to be the case, addressing the other 
terms in Sec. [3j In this state correlations decay exponentially (with inverse correlation 
lengths 7 tA/(/ 3/2) where A is a boundary scaling dimension [12] ) and there is area- 
law entanglement entropy oc clog/3. (We assume that (3 is always larger than any 
microscopic length scale, in units where v = 1.) 

We are therefore faced with the problem of computing correlation functions 


($n(Zn, T n ) ... $2(^2, T 2 )$i{x 1, Ti)) 

in the strip —(3/4 < r < (3/4 with conformal boundary conditions, and analytically 
continuing them to real time t 3 —> it 3 . Fortunately the analytic properties of CFT allow 
this to be carried out. 
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The picture is simpler if we conformally transform the strip: define w = x + it and 
the conformal mapping 

w — y z — i , (6) 

which sends the interior of the strip to the upper half plane HI. A correlation function 
of primary operators <E>j is given, apart from jacobian factors, by a similar correlation 
function in HI which depends on the points Zj and their images Zj in the real axis: 


($ n (z n ,T n ) . • ■$2(x2,T 2 )$l(x l ,T 1 )) (X F({Zj},{Zj}) . 


Moreover, F({zj},{zj}) may be written as a linear combination of conformal blocks 
AkFk({zj}, {zj}), each of which is the analytic continuation to z' 3 = Zj of a function 
Fk({zj}, {z'j}) defined on a branched cover of the whole plane C n . Moreover global 
conformal invariance constrains Fk({zj}, {z'j}) to depend only on the cross-ratios such 
as 

Vn = 0 * - ^){z\ - z , j )/(z i - z')(zj - z'j ), 


with branch points at coincident points where rjij = 0,1, oo. (For the two-point function 
there is only one such cross-ratio.) These blocks can be labelled so that each has a 
well-defined power-law OPE behaviour as some independent set of the r/ij —* o[j3 The 
leading term in the OPE always corresponds to the identity operator, since this limit 
corresponds to the case when all the { zj } and the {zj} are very far from the real axis, 
and each conformal block factorizes into a holomorphic part depending on the {zj} 
and its complex conjugate. This then gives the bulk correlator of the same product of 
operators. 

Something remarkable now happens when we continue to real time Tj —> it 3 . In 
that case 

vjj = Xj — tj , w'j = Xj +tj + if3 /2 , 

and 

z . — ^ j _ _^ e l^/P)(xj+tj) ' 

Note that these all lie on the imaginary axis in the z-plane. However z) is now no longer 
at the image point of z r Suppose now that all the tj are sufficiently large so that 


(xj + Q) - (xji - tj 0 > p, (V(j, j')). 


This corresponds to the intersection of all the past light cones with t — 0 being non¬ 
empty, that is all points have fallen inside a horizon, defined as a future light cone 
originating at some point on t = 0. In that case it may be checked p] that 

mf)), 

f It has recently been pointed out jl5j that for higher-order correlation functions in non-rational CFTs 
the behavior at some of the singular points is not implied by the OPE. This does not affect our later 
conclusion which uses only the 2-point function of twist operators. 
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and we may apply the OPE to deduce that in this limit the leading term is the bulk 
correlation function, that is, the presence of the boundary is irrelevant. If we now reverse 
the conformal mapping, we see that in this limit when all the points have fallen inside 
a horizon, the correlator is equal to that on a cylinder of radius (3, that is, the same 
as at finite temperature. The corrections to this come from higher-order terms in the 
OPE. In the 2 -plane they decay as |z| _A , where A is the bulk scaling dimension of the 
lowest dimension non-trivial operator to which JJ” =1 couples, which has a non-zero 
expectation value near the boundary. This depends on the particular boundary state 
| B). In the w-plane, this corresponds to exponentially decaying corrections of the form 

e -(nA./P)((xj+tj)-(x j '-t j ,)) ' 


2.2. Reduced density matrix. 


The above argument leaves open, for example, the question of other correlators, for 
example those of non-primary operators. A more complete result about thermalization 
would be that the whole reduced density matrix of an interval is close to that of a 
thermal ensemble, once the interval has fallen inside the horizon. This would imply that 
all equal-time correlation functions of local operators with arguments in the interval are 
close to their thermal values. This result we now derive. 

Consider an interval (xi,x 2 ) of length l <C L (in practice it will suffice that 
L/2 — £ 3> j3.) The reduced density matrix of this interval is in general 


Pe(t) 


Tv L _ e e iHt \^o)(^ 0 \e iHt 

Tr L l^oX^ol 


in an obvious notation where Ti'x means a trace over the degrees of freedom in the 
interval X. For the CC state this becomes 


M r ) 


Tr L _* e -H(0O+r)\B)(B\e- H OD-T) 
Tv L e-0/V H \B)(B\e-C/P H ’ 


continued to r = it. 

For real r the numerator of this expression is the partition function on a strip S slit 
along (£, r), the rows and columns of pe being labelled by the fixed values of a complete 
commuting set of local fields above and below the slit. The denominator is the trace 
Tr<» of this, equivalent to sewing up the slit. 

Similarly the reduced density matrix in a thermal ensemble 

_ _Tr L-te~ pH 

p 0 Tvl e ~0H 


is given by the partition function on a cylinder C slit along (£, r), divided by the partition 
function on the full cylinder. Now consider the overlap 

Tr (p e • pp) _ Z(S®C) _ 

((Trp|)(Tr^)) 1/2 " (Z(S © S)Z(C © C)) 1/2 ’ 


(7) 
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Figure 1 . The surface S' ® C in the numerator of 0. The cylinder C and the strip 
S are sewn together along A as shown in the inset. 


where, for example, Z(S ©C) is the partition function on S sewn onto C, such that the 
bottom edge of the slit in S is sewn to the top edge of C and vice versa (see Fig. 1). 
(This construction is similar to that used in formulating the Renyi-2 entropy as a path 
integral.) Note that in the ratio in (J7J) the normalization factors all cancel. 

The surfaces in the three factors in (J7J) are all conifolds: the metric is euclidean 
except at the ends of the interval t where they have a conical points with deficit angle 
47T — 27 t. In a CFT we may view their partition functions as correlators of twist operators 
T evaluated in the product of the CFTs on each component |14j . Thus 

Z<yS®C) = _ {T{xi,t)T{x 2 ,t))s®c _ / g x 

(Z(S © S)Z(C 0 C)) 1/2 «T(X!, t)T(x 2 , t)) S9S {T{x u t)T(x 2 , t)) C0C ) 1/2 ’ 

where now, for example, (T(xi, t)T(x 2 , t))^^ means the twist correlator on a direct 
product of the CFT on S with that on C. 

These twist operators behave in many respects like other semi-local operators in 
a CFT. In particular their correlators enjoy similar analyticity properties. Therefore 
we may take over the earlier arguments applied to correlators of ordinary operators. 
Applying the conformal mapping (JH]) , they are related to correlators of twist operators 
on CFTs on 11 © C, HI © H, C © C respectively. Once the points (xi, r) and (x 2 , r) have 
fallen into a horizon, the cross-ratio p i 2 is 1, and we may apply the OPE. 

The OPE of twist operators has been discussed in [161, lTj. In general it may be 
expressed as a sum of products of local operators in each sheet. In the present case 

T(zi) ■ T(z 2 ) = ^ C kl ,k 2 ^k 1 ^k 2 , 

ki,k 2 

where and are a complete set of operators, arranged in order of increasing 
dimension, in copies of the CFT on (H, C), (H, H), (C, C) respectively, and the C kl ,k 2 
are (calculable) coefficients. The leading term comes from taking = 1, so 
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that, once the ends of the interval have fallen inside the horizon, the right hand side 
of (jHJ) is asymptotically unity. The leading correction comes from the lowest dimension 
non-trivial operators ($&!, < hfc 2 ) which have non-zero expectation values. Since these are 
expected to be primary, and such operators have vanishing expectation value in C, the 
leading corrections are only in the first factor in the denominator of ()8l) , corresponding 
to the most relevant operators which have a non-vanishing expectation values on (H, HI) 
with the particular boundary state | B), as was found above. This gives the main result 
(JTj) of this section. 

Note that if we use the same method to compute the overlap between two thermal 
ensembles with different inverse temperatures f3 and (3', the leading contribution to the 
corresponding ratio of partition functions (coming from the propagation of the ground 
state along the length l for which the two cylinders are joined to give a single cylinder 
of circumference /3 + (3') is 

exp(irc£/6(f3 + (3')) f nc£((3’ — (3) 2 \ 

[exp(7rcf'/6(2 / 3 / )) • exp(7rc£/6(2/l))] 1 ^ 2 V 24 f3/3'((3 + (3')) 

The overlap is therefore exponentially suppressed if j3' ^ (3 and i /3,/3'. This shows, 
by the triangle inequality, that p# cannot be close to a thermal ensemble with inverse 
temperature f3' ^ ft, 

3. Quench from a more general boundary state and the GGE 

In the section we consider a quench in a pure CFT from a more general state given by 
([2]), which we rewrite as 

|^o) tx e~^ /4)HcFT Y[e~^^ k(x)dx \B ), (9) 

k 

where the product is now over all boundary operators excepting the stress tensor. Note 
that we have commuted e~OD) H cFT through all the other terms so it stands on the left. 
Since commutators of H with a given local operator generate others, the effect of this 
is to modify the couplings (3k- Although this expression is motivated by the idea that 
we should include only irrelevant boundary operators with scaling dimension A*, > 1, 
as we argue below it makes sense in some cases also to include relevant ones. 

As it stands, the expression ([9]) is only formal. If expanded in powers of the 
it will involve integrals of boundary correlators of the which will give rise to UV 
divergences if A*, > 1. These may be dealt with in the usual way by first imposing a 
UV cut-off in space, identifying the UV divergent terms by using the OPE, and then 
adding counterterms coupling to the higher dimension operators in this OPE in such 
a way as to cancel the divergences as the cut-off is removed. This has the effect that 
the renormalized f3k appearing in ([9]) then depend on the bare values in a complicated 
non-linear fashion. 

This perturbative expansion should make sense for all irrelevant boundary operators 
for which A*, > 1. For relevant operators one should expect to encounter infrared 
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divergences signaling the crossover to a new boundary fixed point. However in this case 
the finite effective inverse temperature provides a cut-off, and the expansion should still 
make sense as long as (3, that is /3k -C /3 _ ( 1_Afc ). 

Let us first consider the case when the <£*. are all descendants of the identity 
operator, plus their hermitian conjugates. These can be written as linear combinations 
of all possible powers of the components of the stress tensor (T, T) and its derivatives. 
In fact, since the conformally invariant boundary state satisfies 

(T0r)-T0r))|fl)=O, 

there is no need to distinguish between T and T on the boundary. We denote these 
powers and derivatives of T by T^ k \x). They should be normal-ordered as usual by 
point-splitting and removing the singular terms in the OPE. The T^ have the property 
that their boundary scaling dimensions A*, are integers > 2. In fact, they should be 
even integers for the state in (j9|) to be invariant under reflection x —>■ —x. Note that 
there is in general more than one such descendant operator with dimension A*,. 

As remarked above, there is no need to distinguish between T( k \x) and x ) on 

the boundary. However they may each be viewed as the boundary values of holomorphic 
and antiholomorphic operators T^ k \w) and T^ k \w) which are different in the bulk. 
In the path integral in the euclidean slab geometry, each exponential factor in (l9j) 
corresponds to inserting a term 

/3 k J T (k \x)dx = \p k J T i - k \w)dw + \p k J T [k \w)dw 

along contours just above and below the lower and upper boundaries r = ±/3/4. There 
is a contour for each k, and they should be ordered in increasing k. The main point 
now is that, because of the (anti-)holomorphicity, each of these contours may be freely 
deformed into the bulk, as long as they do not cross each other or the arguments of 
local observables in correlation functions. This corresponds to the statement that the 
charges which are the space integral of the T^ commute with the hamiltonian. On 
transforming to the z-plane, the insertion becomes (for each k) 

±(3 k (-i(3/i r)) 1 -^ J T^ k \z)z Ak ~ 1 dz + • • • + c.c., 

summed over two contours: one from z = 0 to positive real +oo, the other from z — 0 
to real — oo. The omitted terms in the above correspond to the contributions of other 
(holomorphic) descendants, which appear because the T^ are not primary. However 
these will disappear when we reverse the mapping going to the cylinder. 

We now continue the arguments of the local observables to real time. As before, 
we argue that at late times the Zj move off to — zoo so that effectively the boundary 
disappears. However, the boundary perturbations given by the above contour integrals 
remain. In the absence of the boundary we can now rotate one contour into the other, 
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for example choosing them both to run from z — 0 to —ioo. This gives rise to a relative 
factor 

( e -W2)A fe + ( e W2)A fc = 2 cos ( 7rAfc /2). (10) 

These manipulations are illustrated in Fig. 2. 

Although Afc is, so far, an integer, we have kept the more general expression for 
later use. Note that only even integer dimensional charges contribute: the odd ones are 
odd under parity and are zero acting on the initial state. 

On transforming back to the w-plane (which is now a cylinder) we find an insertion 
in the path integral action 

A,cos(7rA fc /2 )H ( - k) , 

k 

where 

/**> = /(**>(,, T)+ 

along r = f3/ 2. However, because commutes with the hamiltonian, this contour 
could be along any constant imaginary time line which does not separate any of the 
arguments Wj of the observables. In addition, the contours for different k should be 
correctly ordered, reflecting the fact that the do not in general commute among 
themselves. 

Thus the correlation functions at times after all the points have fallen within the 
same horizon are given by a path integral with weight 

e ~8cFT~J2k Acos(vrA fe /2)H( fc ) 

where the first term is the usual action for the CFT. This is just the path integral 
formulation of the (non-Abelian) GGE. Note that odd-dimensional charges do not 
contribute. 


3.1. Physical consequences of the GGE. 

In what ways does the post-quench GGE differ from the Gibbs ensemble in measurable 
quantities, for example the correlation functions of local observables? Recall that in a 
Gibbs ensemble in a CFT, the stationary 2-point function of any scaling operator decays 
as 

($(xi, t)$(x 2 , t))y ~ e - 2 ’rA *|* 1 -* 2 |/P ? 

for \x\ — xf\ 3> (3, where A$ is the overall scaling dimensions of <f> (assumed to be a 
scalar). The reason for this is that the generator of translations along the cylinder is 


h = (1/2tt) / T xx dr = (1/2tt) / (T + T)dr = (2tt//3)(Z 0 + l 0 - (c/12)), 


where 

T(t) = ) 2 ^(i, - , T(t) 


(2 


27r iqr / (3 


( 12 ) 
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Figure 2. The contour manipulations leading to m- (a) initially the perturbing 
operators are integrated along the dashed lines adjacent to the upper and lower 
boundaries of the strip; (b) using (anti-)holomorphicity they are moved into the bulk; 
(c) the strip is mapped to the upper half-plane: the correlator of local operators 
(indicated by black circles) depends on those points and their images; (d) after 
continuation to real time and evolving until both operators fall within the horizon, the 
image points move off to infinity: at this point the boundary (now indicated by a dotted 
line) effectively disappears and the correlator is the same as in the full plane; (e) the 
two parts of the dashed contour may now be rotated together and their contributions 
(which differ by a phase) combined; (f) on reversing the conformal mapping, the full 
plane becomes a cylinder with the insertion of a defect line along which the conserved 
currents are integrated. 
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(we use lower-case letters to distinguish these from the hamiltonian H and the associated 
Virasoro generators L q , L q .) The eigenvalues of Q + lo are the total scaling dimensions 
A(j>. 

In the GGE, the additional terms contribute to the action 

c ^ 

V'/3 fc cos(7rA fc /2) / (T (fc) (x, r k ) + T (fc) (z, r k ))dx . 

k Jo 

where the Tk are ordered in increasing k. 

In terms of the generator of translations along the cylinder, there is now an 
additional term 

h ->■ h+6h = /1+/3' 1 ^/3 fc cos(7rA fc /2) f (:T (fc) (r fc +r): + : r (A:) (r fc +r) :)dr .(13) 

k Jo 

We have averaged over the overall center of mass of the {t*,} to emphasize the 

stationarity of the ensemble. This expression may be written in terms of the L,L. 

y y _2 

As an example, the lowest (dimension 4) operators correspond to : T 2 : and : T :. 
Inserting the Fourier decomposition (I12p . we have, before normal ordering 

5h = /5 2 (2tt//3) 4 J^( l q - ^dqo )(Z_, - ±5 - q0 ) + c.c. 

<? 

In this case, the effect of normal ordering is to move all the l q with q > 0 to the right. 
Subtracting off the singular terms in the OPE corresponds to neglecting the commutator 
when we do this. Thus 

Sh = (3 2 (2 i t//3) 4 ( 2 ^ l_ q l q + (Z 0 - ^) 2 ) + c.c. 

V q >0 ) 

The terms with q > 0 annihilate any highest weight state, corresponding to the ground 
state of h or the insertion of any primary operator at infinity. Thus the partition function 
of the GGE is now 


Z ~ exp L ((nc/6/3) — 2/? 2 (27r/ / d) 4 (c/24) 2 ) . 

Thus the energy density — (1/L)dp logZ is modified, and it is no longer simply 
proportional to the entropy density (1 /L)(/3dp — 1) log Z. 

The 2-point function of a scalar primary operator now decays with a correlation 
length 

($ %ge r 1 = (27rA*//3) + /3 2 (27r//3) 4 ((A* - ^) 2 - (^) 2 ). (14) 

Thus in the GGE, the effective temperature as measured thermodynamically, or by a 
correlation function, and therefore in terms of fluctuations, depends on the observable. 

Higher dimension operators may be dealt with in a similar way, although the algebra 
becomes more complicated. The form of the next few normally ordered operators has 
been given in Ref. [8j. The general picture does not change: the correlation length in 
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the GGE will be expressed as a power series in the scaling dimensions of the particular 
observable. 

Note that this will also have an effect on dynamical quantities. Since the time 
evolution remains governed by Hcft, there will still be a sharp horizon (rounded on 
the scale of (3), and the spatial decay of a 2-point function in the stationary state for 
\x\ — x 2 | < 2 vt should match across the horizon with the time decay of the product 
of the 1-point functions for \x\ — a^l > 2 vt. Thus a 1-point function (which does not 
vanish in the initial state) should decay as 

(®(x,t)) ~ e ~ vt/ ^GE . ( 15 ) 

This has been verified for the Ising model [18] with an explicit form for S,gge- 

3.2. Perturbative derivation. 

The above argument relies on distorting space-time contours in the exponential 
weighting of the path integral. Since this is somewhat heuristic, it is therefore reassuring 
to see that similar consequences can be seen within straightforward perturbation theory 
in the parameters (3k- 

We illustrate this by considering the modified decay of the above 1-point function 
©, and the simplest example of a T 2 operator perturbing the initial state. This may 
be evaluated by continuing to imaginary r the sum of integrated correlators 

(1/n!) 

n =0 

where the integrals are along the upper and lower edges of the strip. (As discussed above, 
these must be regulated with a UV cut-off and the divergent terms subtracted.) The 
correlations functions are in the CFT with the unperturbed boundary state. As such, 
they may be estimated by transforming to the z-plane El as before. The correlator 
depends on (0, r) and its image, which as before maps onto the imaginary z-axis, 
although not at conjugate points once the continuation to real time is made. The 
points Uj map onto points Xj = ±_^ 2 ' K IO u i on the real axis. 

Once again, points separated by distances > OQ3) in the w-plane will become 
exponentially ordered in the z-plane, and, apart from these cross-over regions, we may 
apply the OPE. 

Let us consider the case n — 1 as an example. This involves the 2-point function 
of the bulk operator <f> with the boundary operator T 2 . Since the form of this follows 
from global conformal invariance, at this order we may initially consider a more general 
boundary operator $ 2 , with scaling dimension A 2 , rather than T 2 . 

The first order correction is 



dui... du n ($( 0, r) T 2 (ui, ±|)... T 2 (u n , ±^)}du 1 ...du r 
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As before we conformally map this into the upper half plane by z — ie 2nw ^, giving 



\dw/dz\ 2A ‘ s> \du/dx\ 1 


dx 

( i(z — z)) 2A *~ A2 (z — x) Az (z 


x) A2 ’ 


plus a similar contribution from the other boundary, which means replacing x —> —x. 
We have used the explicit form of the 2-point function in H (which has the same form 
as the 3-point function in full plane.) Here & $< g 2 is a bulk-boundary OPE coefficient. 
Continuing t —^ it as before, we have ~ (/3/27 t), |du/c7r| rs j (/3/2nx), giving 


-^2^5 2 (2tt/^) 


2A$—A 2 


( e 2 W/3) 


A 2 —2A«|) 


d 2ttu//3 


2ttu/(3 _ ^g—2'Kt/(3^j^g27Tu/(3 _j_ ^g27rt//5 ^ 


A 2 


du , 


plus the complex conjugate contribution from the other boundary. 

For u < —t — 0(/3) we may approximate the integral giving 

-^^|, 2 (27r//3) 2A *- A2 (e 27rf//3 ) A2 - 2A$ f (e 2 ™/P) A2 du = —/3 2 6 <I>3>2 (2tt//3) A *” A2 

J — OO 

We get the same contribution from u > t, and from the other boundary for \u\ > t. 
These all give an O(fo) correction to the pre-factor of the 1-point function. 

However, from —t<u<t we have 

-fobs $2 (27r//3) 2As ” A2 (e 2 ^ //3 ) A2 - 2As [ 0 ie 2nt/2fi )~ A2 du. 


Adding the contribution from the other boundary we get 


-fobs $ 2 cos(7tA 2 /2) t ( 27r / / g) 2A *- A2 (e 2 ^)- 2A * , 

so, relative to the O^fo) term, we get a factor 

1 - fobsl 2 cos(7rA2/2)(27r//3 ) _A2 1 . (16) 

Therefore the first-order correction appears to give a larger contribution as t —» 00 
than the zeroth-order term. The reason that it is oc t is physically clear from the 
quasiparticle picture: at time t the perturbation of the initial state is felt only from 
regions with |u| < t, by causality. Within this region, however, the integrand is 
approximately constant. Therefore the contribution is oc t. For the same reason, 
we then expect the leading behavior of higher-order contributions to behave like 
(—fo) n t n /n\, so that the first-order correction simply exponentiates. This can be seen 
more clearly in the z-plane: the leading contribution now comes from the region where 
\z\ <C X\ <C X 2 <C • • • <C x n <C \z\. In this limit we may use the bulk-boundary OPE 
repeatedly, each time getting a factor Note that the factor cos(7tA 2 /2) in ([3]) also 
naturally arises in this approach. 

Specializing to the case $2 = T 2 , we have by the conformal Ward identity that 
bs t T 2 — A|, which then confirms the term 0(A|) in (fl4l) . (The other term oc c arises 
from normal-ordering.) 
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The result (fIBl) and its exponentiation have a simple interpretation in the 
quasiparticle picture: the perturbing operator acting on the initial state gives rise to 
additional quasiparticle pairs. Only those coming from the past light cone of (x, t ) 
affect {$(x,t)), and those coming from initial points (a;, 0) separated by distances > f3 
are incoherent, so their effects exponentiate. 

3.3. Non-integer exponents and semi-local conservation laws 

The above argument assumed that the operators acting at the boundary have 

(even) integer dimension. However, in most CFTs and for most conformal boundary 
states there will also exist operators with non-integer dimensions. One may ask whether 
these also lead to a GGE in the stationary state. For the above argument to apply, we 
have to be able to think of these boundary operators in some sense as the limits to the 
boundary of bulk holomorphic and anti-holomorphic operators $fc(z) and $*,( 2 ). We 
argue that, at least for rational CFTs, this is indeed the case. Recall that in this case 
all primary operators (under Virasoro or some extended algebra) have descendant null 
states which implies that their boundary correlators (^(aq)^^)...) satisfy systems of 
linear differential equations with singular coefficients at coincident points. The solutions 
then have power law singularities at these points, but may otherwise be continued 
uniquely into the upper (or lower) half-plane as holomorphic (or anti-holomorphic) 
functions. These may be viewed as the (multiple-valued) correlators of holomorphic (or 
anti-holomorphic) operators. If the boundary operator has boundary scaling dimension 
A*,, these bulk operators have dimensions (Aj.,0) and (0, A*,) respectively. 

Another way to construct these bulk operators, which works in principle for any 
CFT, is to cut out a small disc of radius e around a bulk point (z,z), then insert <!>*. 
at the point z + ee l6 on the boundary of the disc. Correlators of the bulk holomorphic 
operator are then defined as weighted integrals of those of the boundary operator: 



(M*) ■■■) = lim / e~ lAk0 ($ k {z + ee ie ).. )d0 . 


This resembles the construction of discretely holomorphic bulk operators in lattice loop 
models [19]. 

Having defined these bulk operators, the previous construction of the GGE follows, 
at least formally. However the associated charges Qk = J ^k(u,r)du + J &k(u,r)du 
do not quite commute with the hamiltonian He ft- This is because if we consider a 
closed contour which is the boundary of a long rectangle {—\L < u < |L, T\ < r < r 2 ), 
the contributions from the end pieces do not cancel for periodic boundary conditions 
since they carry different phases e ±lAk . They would in fact cancel for suitable twisted 
boundary conditions on the fields of the theory. Thus we may think of the action of a 
single Qk as switching the boundary conditions between different sectors of the theory. 
It is a parafermionic charge. If A k is rational, with lowest denominator M, it is only 
after acting M times that we return to the original boundary conditions. Thus, in the 
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GGE expression 

Tr e~^ HcFT e~^ k cos C A k/ 2 )Qk 

k 

the trace is supposed to project onto only those terms in the expansion of the exponential 
containing integer powers of /3jf. 

3.3.1. Example: the transverse Ising model in a longitudinal field. As an example, we 
consider the CFT description of a quench in the Ising chain with hamiltonian 

-ffIsing = - d ^ (7,J a ] — h x ^ (Tj — h z ^2 a j > 
j 3 3 

from the ground state with h x ^> J and small longitudinal field h z -C J to the critical 
point h x = J, h z = 0. This has not been considered in the literature because the model 
is no longer solvable in terms of free fermions when h z 0. The appropriate conformal 
boundary state when h x J and h z = 0 is that corresponding to free boundary 
conditions on the crj. This state supports one primary operator of scaling dimension 
S — \ which is interpreted as the scaling limit of the local magnetization cr z . Turning on 
a small h z is equivalent to perturbing the boundary state as described above. Although 
in this case the perturbation is relevant, as was argued above, we may still consider it 
as small as long as h z <C /W 1 / 2 ~ (h x /J — l) 1 / 2 . Note that the expansion in terms of 
correlators of o z on the boundary contains only even powers of h z because of the Z 2 
symmetry of the state. 

The holomorphic and anti-holomorphic extensions of this boundary operator are in 
this case well understood - they are nothing but the fermions (i/j( z), ijj(z)) of the bulk 
Ising CFT. Thus the GGE in the case should contain fermonic charges (Q, Q) which, as is 
well known, act to switch between periodic and anti-periodic boundary conditions on the 
Ising spins a z . If we consider a fixed choice of boundary conditions, for example periodic, 
only even powers of these charges enter the computation of correlation functions of local 
spins in this GGE. This is completely consistent with the fact that only even powers of 
h z enter the perturbative expansion. 

4. Perturbations of the CFT 

In this section we consider what happens in the more realistic case when the time 
evolution is governed by a hamiltonian H which differs from that of a pure CFT by 
the addition of couplings to irrelevant operators. While these give only corrections 
to scaling at the zero-temperature critical point, their precise influence on the quench 
dynamics has not so far been investigated. Studies in integrable models suggest that 
the quasi-particle picture remains valid, with due allowance for the different dispersion 
relation. 

While it is possible to study these effects perturbatively in the couplings to these 
irrelevant operators, similarly to the case of the deformed initial state in Sec. |3l this 
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similarly breaks down at late times, and, unlike the former case, it is not obvious to see 
how to resum the series. Instead we take a different approach, which, as we will show, 
agrees with perturbation theory at low orders. Unfortunately this approach only works 
for irrelevant operators which are descendants of the identity, that is, polynomials in 
the components (T, T) of the stress tensor and their derivatives. In fact, most of the 
analysis will be directed at the dimension 4 operators T 2 , T and TT, arguing later 
that the general picture should remain valid for higher dimensional operators if weakly 
coupled. However the dimension 4 operators are important, as they are generated in 
almost any quantum critical hamiltonian with a linear low-energy dispersion relation 
by integrating out higher energy degrees of freedom up to the scale of the UV cut-off J§] 
The first two terms break Lorentz invariance, and would be generated, for example, in 
a lattice model. However, as discussed in the previous section, their space integrals in 
fact commute with the hamiltonian and they preserve the integrability of the model. 
However the third term, while preserving relativistic invariance, introduces right-left 
scattering and does not commute with This may be seen by quantizing the 

CFT on a circle of length L , whence 


f>L _ pL _ r L 

/ T 2 (x)dx oc / T~(x)dx LkL_k, / T(x)T(x)dx oc LkL & 

j 0 7 „ J 0 7 „ J 0 1- 


The first two commute with Hcft oc L 0 + Lq, while the last does not. 
Let us therefore consider the perturbed hamiltonian 


H 


H, 


CFT 


hx 


J (T 2 + T 2 )dx — n 


TTdx , 


which corresponds to a euclidean action 


S E 


qE 1 \ 

°CFT 2 A 


J (T 2 + T 2 )dxdr - /a 


TTdxdr. 


Note that we do not normal order, but will shortly introduce an explicit UV regulator. 
The coupling constants (A,/i) have the dimensions (length) 2 . The additional terms, 
quadratic in T and T, may be generated by the usual trick of introducing auxiliary 
complex conjugate fields (£,£), with spins ±2, in the functional integral: 

S E = S e ft + J (J 1 & - |A(e 2 + e))dxdr + J (£T + QT)dxdr (17) 

In Feynman diagram terms the interaction may be thought of as proceeding by the 
exchange of (£,£) ‘particles’ which couple linearly to T and T. A UV regulator may 
be introduced by adding a kinetic term oc f (dS)(dS)dxdT. However, this makes sense 
as a unitary theory in Minkowski space only if the eigenvalues of the (mass) 2 matrix 
(fi ± A) -1 are positive. One might assert the freedom to choose the signs differently, 

§ In a free theory, with has conserved 17(1) currents (J, J), the curvature of the dispersion relation also 
results in terms like dJdJ. 

|| However it may still be integrable, see for example BU . 
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but these would not result from the integration over unitary relativistically invariant 
higher energy fluctuations, at least of a scalar nature. As we shall see, the choice of 
the ‘wrong’ signs would lead to superluminal speeds of propagation in the low-energy 
effective theory. This, of course, is not a problem if this theory is describing a lattice 
model which does not have to possess a relativistic ultraviolet completion. 

Although we shall continue to work in euclidean space for the time being, continuing 
to real time only at the end, it is interesting to note that in Minkowski space cm 
corresponds to 

S M = S% FT + \\ j (T 2 + + T 2 _)dxdt + /i J T ++ T__dxdt , 

where we have introduced light-cone coordinates x ± = x±f. It is instructive to consider 
the case of a free scalar held 0, for which T±± oc ( d±cj )) 2 . In momentum space the in¬ 
teraction vertices are then — iX(ki + k 2+ k 3+ k 4+ ), — i A(/ci_/c 2 -& 3 -& 4 -), —ifi(ki + k 2 +k 3 -k 4 _). 
H ± A > 0 then corresponds to an attractive interaction. (This may also be seen by com¬ 
paring (H3 with case of a simple non-derivative interaction S = Sf ree — /a J cf) A dxdT. 
The ‘usual’ sign with /a < 0 is repulsive.) 

Note that the couplings A, /i have dimension (length) 2 . In what follows we assume 
that they are small compared the dominant length scale of the problem which is /3 2 . 
This means that we can assume that the (dimensionless) fields (£, £) are -C 1. 

The main point now is that the last term in (1771) may be reinterpreted, using the 
definition of the stress tensor, as the response of the action to a small change in the 
metric 5g = 5g ww — £, 6g — 5g ww = £, which themselves may be viewed as gaussian 
random fields with covariance 

E[hg(w, T)dg(u', r')] = A 5(u — u ', r — t') , 

E [5g(u, r)5g(u ', r')\ = A S(u — u' : r — t') , 

E [6g(u, r)5g(u t')) = g,5{u — u', r — r'). 

If we introduce a UV regulator as above, the delta functions are smeared over a length 
scale 0(A 1 / 2 ,/i 1 / 2 ). 

The deformed metric is 


ds 2 = dwdw + 5gdw 2 + 5gdw 2 . 

To lowest order (only) this is equivalent to a coordinate change ds 2 = dwdw where 

dw = dw + 5gdw, dw = dw + 5gdw . 

Again to lowest order, this implies that the correlation function of any product of local 
observables evaluated with the perturbed action (11711 is given by the expectation value 
over the random fields of the same product in the pure CFT at shifted values of the 
arguments: 

(0( w j,Wj))s CFT +SS ^j)) Soft] ■ 
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where E[-] means the average over the gaussian fields. 

We now apply these results to the quench scenario. Recall that a point (. Xj,tj ) in 
Minkowski space is mapped to 

Wj = Xj — tj , Wj = w'j = Xj + tj + i/3/2 , 

on the cylinder. Thus the effect of the coordinate shift in light-cone coordinates is 

dx + = dx + + Sg dx- , dx- = dx- + Sg dx + . 

Note that we have avoided going into Minkowski signature until this last step: although 
a similar argument should be possibly directly in real time (thereby being applicable 
to the effect of the perturbation on more general time evolution besides that of quench 
scenario) in practice we have been unable to resolve tricky questions of reality in the 
gaussian measure. 

We now discuss some simple consequences of this picture. 

4-1. Renormalization of propagation speed. 

In the shifted (tilded) coordinates the correlations are to be evaluated in the pure CFT, 
so signals propagate with speed v (which so far we have taken to be 1.) The left- and 
right-moving null geodesics are given by dx + = 0 , = 0 respectively, that is in the 

original coordinates 

dx + = — Sg dx _ , = —Sgdx + . (18) 

The first effect is due to the fact that, in the quench scenario, Sg and Sg have non¬ 
zero expectation values. This is because (T) and (T) are both non-vanishing, see (J5]) . 
From 03 we that see that at the saddle-point (Sg) = (£) = —(p + A)(T) 7 ^ 0, and 
similarly for the complex conjugate. To this approximation (neglecting fluctuations), 
the null geodesics become dt = ±(1 + 2(8g))dx, so the average speed of propagation is 
renormalized. Since (T) = (T) < 0 in the initial state (and thereafter), (Sg) and (Sg) 
are both positive if p + A > 0, and so the speed is reduced. 

More directly, we see from (TT7T) that there are terms in the hamiltonian 

SH =- I (A (T)T + A (T)T + p(T)T + p(T)T)dx . 

Recalling that (T) = (T) = — 7 rc/ 12 / 3 2 , we see that this term is 

SH = (A + p)('kc/Q(3 2 )H . 

This is equivalent to a rescaling of the time coordinate t —y t( 1 + (A + p)(nc/6/3 2 )). On 
the other hand, a similar argument applied to the generator h of translations along the 
cylinder or strip yields the opposite sign for the shift (because T xx = — T rr in the CFT) 
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which is equivalent to a rescaling x —* x(l — (A + /r)(7rc/6/3 2 )). The propagation velocity 
is therefore renormalized according to 

v —> v(l — 2(A + /t)(7tc/6/S 2 )) . 

If A + /i > 0 this corresponds to a reduction, which is physically reasonable: the 
attractive interaction causes a positive time delay as the quasiparticles scatter from 
the rest of those in the finite density sea. Note that this attenuation is greater at 
higher effective temperatures as expected. If we had chosen the other signs for the 
couplings (corresponding to a non-unitary interaction with the higher energy scales) 
this would have led to an increase in propagation speed with increasing temperature. 
In a microscopically relativistic theory this would be impossible. Similar observations 
(of couplings with the wrong sign implying superluminal propagation) have been made 
in the context of unification in particle physics m- 


4-2. Diffusive broadening of the horizon 

We now return to ([TED an d include the short-range correlated random nature of (5g, 8g). 
For convenience we assume that the rescaling of x and t implied by the fact that they 
have non-zero expectation values has already been performed. If we integrate (fT8l) to 
lowest order in the couplings, the equations for a right-moving null geodesic starting at 
(x + = 0,x_ = £_(0)) is 

j-x+ 

x_(a; + ) = x_(0) — / 8g(x' + ,x_(0))dx , + . 

Jo 

If the integrand were delta-correlated only in x ' + , this would lead to Brownian motion 
for x_(x + ). However in fact 

(8g(x' + , x'_)8g(x+, x"_)) = A5(x' + — x" + )8{x'_ — x "_), 


so that the effective diffusion constant will be divergent. Even if a UV regulator is used 
to modify the correlator of the random field, the result will depend strongly on this. 
Thus individual null geodesics will behave rather wildly on the scale of the UV cut-off. 

However, recall that the horizons in the pure CFT have a width 0(/3). If we define 
the mean coordinate of the horizon by the random variable 


rhi 2 

X_(x + ) = /3 _1 / x_(a; + )(ia;_(0), 

J-/3/2 


this has variance 

<X_(a; + ) 2 > = ( \/(3)x + . 


(If we use a more realistic profile for the unperturbed horizon, the coefficient will change 
but not the functional dependence.) Thus we see that the effect of the T 2 and 
operators is to diffusively broaden the horizon to a width 0((\t/fJ) 1 ^ 2 ). However this is 
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still smaller than t so on this scale the horizons remain well-defined. Note also that the 
TT term, which describes left-right scattering, does not contribute to this effect. This 
can be seen in the example of the free boson discussed above. Since the scattering is 
proportional to (k\+k 2 +kz-k±J), it vanishes on both left- and right-moving light-cones 
and so does not modify them. This is consistent with a pure TT interaction preserving 
exact Lorentz invariance. 

4-3. Correlation functions. 

We now discuss how the coupling to the random metric affects the asymptotics of 
correlation functions. First consider an equal-time correlator (^(aq, t)$(x 2 , t)) in the 
stationary state, that is when 2 1— |aq — xf\ 3> /3. In the unperturbed theory its decay is 
given by (£ / f) _1 = 2ttA^//5. Consider the space-like geodesic between these two points. 
Along this we have 

dx = \{dx + + ete_) = dx + 5gdx_ + 5gdx + ~ dx + ( 5g + Sg)dx , 
since to this order dt can be neglected. Thus 

mx u t)*(x 2 ,t)) H+ 6H ~ E [ e - 2 *A*l*-«l//»] , 

where the expectation on the right hand side is with respect to the gaussian random 
fields. From the above we see that X\ — x-i is a Brownian motion, but, once again, 
because the fields 8g, 5g are defined at equal times, the effective diffusion constant 
depends strongly on the UV cut-off. Again the resolution is not to consider the x- 
cooordinate at a particular imaginary time r, but averaged over the cylinder. This is 
equivalent to computing the stationary behavior. This variable has a finite variance 
2(A + /a)/(3, and so we find 

($(xi ,t)$(x 2 ,t)) H +5H ~ exp j | Xl - ^ . 

This result agrees with first-order perturbation theory, since (A$|T 2 ,T ,TT |A$) = 
(2vr//3) 4 A|. 

4 . 4 . Higher order effects. 

Finally we consider the effects of taking into account higher orders in A and /i, as well 
as higher order descendants of the stress tensor in the quench hamiltonian. The 
latter terms can still be described by a coupling to a random metric with non-gaussian 
correlations, which should not change the overall picture if they are small. 

In general, since the metric is 2-dimensional, it is locally conformally flat, so it is 
always possible to choose a local system of coordinates x± so that the metric has the 
form 


ds 2 = dx + dx_ . 
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However, in higher orders this will correspond to non-zero curvature R oc <9 + <9 _ 0. 
This first arises at second order, when R oc d + d~{5g + 5gD). Thus it first appears as 
a consequence of (potentially non-integrable) L-R scattering. In fact, pure LL and RR 
scattering is always equivalent to a simple coordinate change x± —> x± and introduces 
no curvature. Thus in the latter case we expect to see only broadening of the horizon, 
possibly non-gaussian. On the other hand, strong negative curvature effects could lead 
to other effects such as the chaotic divergence of geodesics. However this seems difficult 
to analyze quantitatively. 

Note that while this discussion is reminiscent of 2d quantum gravity, the physics is 
very different - the ’graviton’ here has a heavy mass 0((p± A) -1 / 2 ). In fact, the physics 
is more similar to the scattering of light by gravitationally bound dark matter in the 
distant universe. 

5. Discussion. 

In this paper we have considered various extensions of the initial results in CC [1J 
concerning the quench to a critical point in a 1+1-dimensional CFT from a state with 
short-range correlations. First, we showed that once a finite interval falls inside the 
horizon its reduced density matrix becomes exponentially close in L 2 norm to that of 
a thermal state, generalizing the results of CC which applied only to correlators of 
primary operators. It would be straightforward to extend this to other L p norms using 
the same techniques. Second, we showed that deformations of the CC initial state by 
action of integrals of local boundary operators should lead to a (non-Abelian) GGE with 
corresponding local conserved charges. However, it is in general necessary also to include 
parafermionic charges. Finally, we showed that irrelevant perturbations of the CFT 
hamiltonian involving powers of the stress tensor lead to dependence of the propagation 
speed on the initial state, and to a diffusive-like broadening of the horizon. Although we 
have treated these three topics independently, they could clearly be combined, at some 
technical expense, into a unified picture, considering, for example, the convergence of 
the reduced density matrix to that of a GGE in the case of a general initial state, in a 
perturbed CFT. 

There has been some discussion in the literature as to exactly which conserved 
charges should be taken into account in the GGE ED The results from lattice studies 
initially suggested that only ‘local’ operators should be included, defined as those which 
involve lattice spins over a finite range of lattice spacings. However these have been found 
to be incomplete in some cases, and this appears to be associated with the existence 
of bound states [23]. However, recently this problem has been repaired for the case of 
Heisenberg chains by the addition of ‘quasi-local’ charges |22j . It has been similarly 
argued within integrable field theory that the usual local conserved charges are not 
sufficient to describe the GGE [2HJ . 

The case of a rational CFT, however, appears to be different: in a finite but large 
system the eigenstates of the hamiltonian are in correspondence with the states of a 
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finite number of Virasoro modules (or those of an extended symmetry) and these states 
are in bijection with the Virasoro conserved charges (i.e. those made up of the stress 
tensor and its derivatives.) Thus, giving the expectation value of the conserved charges 
in any initial state should specify that state uniquely. (In fact, because of the existence 
of null states, the Virasoro charges are over-complete.) This argument suggests that 
for a rational CFT, the non-Abclian GGE is sufficient. One might ask whether less 
information is needed to specify the reduced density matrix pi in a subsystem of length 
i <^L. In fact, by considering propagation along the cylinder, it can be seen that pe is 
determined in this limit by only the ground state of the hamiltonian h + Sh in (TT3|) . not 
its full spectrum. Since this hamiltonian is defined on a system of length /3 <C Q one 
might argue that this requires less information that specifying the initial state. However, 
this argument compares two infinities and it seems difficult to make more precise. 

Nevertheless the values of the charges in the commuting sub-algebra given in |8], 
are insufficient to specify pe, since it retains perfect memory of all the left- and right- 
movers in the initial state in regions of length l. This feature appears to be related to 
the exactly linear dispersion relation in a CFT. A similar result has recently been found 
by Sotiriadis [9] for the case of a free boson evolving from a non-gaussian initial state. 
Although this is a CFT, it is peculiar in possessing U(l) conserved currents (J, J), out 
of which local conserved quantities may be formed by taking products and derivatives. 
(Alternatively, the conserved quantities are the (non-local) number operators of the 
momentum modes of the held. The linear dispersion relation appears to play a crucial 
role for this model, since it was shown ra that for a massive free boson memory of the 
initial state is lost and the number-operator GGE is sufficient. 

However, our argument does not necessarily extend to CFTs where the central 
charge is greater than the number of primary conserved currents, with an infinity of 
primary holds. These CFTs are believed to be non-integrable, or chaotic, in a well- 
defined sense ra- Nor does it apply to massive integrable held theories viewed as 
perturbed CFTs, since in this case only some of the Virasoro charges remain conserved. 

The parafermionic operators we claim that one must also introduce are better 
described in CFT language as ‘semi-local’: they depend on single space-time coordinate 
but depend non-locally on the state of the system described in the basis of local 
operators. Nevertheless they appear in the spectrum of the CFT, for example on a 
cylinder, if suitable boundary conditions are imposed, and in some cases, such as the 
Ising fermion, are more fundamental than the local helds of the theory in the sense 
that the latter may be built out of them. It is interesting that it should be in principle 
possible to access their correlation functions in the quench scenario according to the 
above arguments. 

It would be interesting to extend our analysis of the non-integrable TT perturbation 
of the CFT to other situations, for example an inhomogeneous quench [25j, where there 
are non-zero energy currents, or to the steady state currents set up by a non-zero 
temperature gradient EH- Recently this problem has been addressed from a different 
point of view [28]. 
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